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Abstract. Recent experiments with ultra-cold atoms have demonstrated the 
possibility of realizing experimentally fermionic superfluids with imbalanced spin 
populations. We discuss how these developments have shed a new light on 
a half-century old open problem in condensed matter physics, and raised new 
interrogations of their own. 



1. Introduction 

Since its discovery nearly a century ago by Kammerling-Onnes, superconductivity has 
remained one of the most active field of research in physics. The study of this dramatic 
feature of solids at low temperature has been the source of countless new concepts 
and applications, from the discovery of its microscopic origin by Bardeen, Cooper 
and Schrieffer (BCS), to the development of MRI in medical imaging or the invention 
of SQUIDs (Superconducting Quantum Interference Devices) which constitute the 
most precise magnetic field probes. Nowadays, they are used in the development of 
quantum computing for the experimental realization of q-bits, and the understanding 
of the microscopic origin of high critical temperature superconductivity remains one 
of the most famous open problem in physics. 

Recently, with the observation of superfiuidity in ultra-cold Fermi g ases, a 
promising connexion has been made between atomic and solid state physics ^ 0. 
In particular, the experimental exploration of the phase diagram of atomic fermionic 
superfiuids as a function of the strength of attractive interactions has confirmed the 
existence of a smooth crossover between the BCS regime of weakly attractive fermions 
and the Bose-Einstein condensation (BEC) of deeply bound pairs proposed in the early 
80's by Leggett [3j, Nozieres and Schmitt-Rink [T. Another intriguing issue recently 
addressed experimentally is the fate of a superconductor when the spin populations 
are imbalanced. Indeed, in the standard BCS theory, superconductivity arises from 
Cooper pairing of opposite spin fermions, and is therefore sensitive to a mismatch 
between the Fermi surfaces of the two spin species. This fundamental question has in 
fact been raised in various branches of physics, including exotic superconductivity [5] 
in heavy fermions, organic compounds or cuprates, and nuclear physics [6 , each field 
bringing its own pecuHarities to the topic. 

Mean field theoretical foundations on this issue were laid in the 60 's but the phase 
diagram of a polarized Fermi gas remained unexplored experimentally until recent 
work in Rice University, MIT and ENS on ultra-cold Fermi gases. In what follows, we 
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present a brief account of the experimental findings made by these groups. We wih 
show how they triggered new theoretical ideas and led to a more refined understanding 
of these systems. Note that this review voluntarily focuses on aspects closely related to 
experiments on ultra-cold Fermi atoms and leaves out more speculative issues. More 
information can be found in related reviews, in particular [H [2j [71 [8] • 

2. The founding fathers 

A significant progress towards the understanding of superconductivity was achieved 
when Cooper [9j realized that an arbitrarily weak interaction pairs electrons with 
opposite spins in the presence of a Fermi sea. Bardeen, Cooper and Schrieffer then 
proposed a variational form for the ground state wavefunction consisting of a coherent 
superposition of these Cooper pairs [10 . This wavefunction lies at the core of the 
standard BCS theory and takes the form 

k 

where |0) is the vacuum and and are variational parameters satisfying the 
constraint uf^-\-v'^ = 1. c^^c^^^ creates a Cooper pair with zero total momentum and 
the ratio Vk/uk is the wavefunction of a single Cooper pair in Fourier space (Fig. [TJa). 
In the standard BCS treatment, the self-consistency of the theory is applied to the 
order parameter A which sets the gap for single-particle excitations. The complete 
pairing of spin t electrons with spin I electrons in the BCS wavefunction ([T]) requires 
the spin populations to be exactly balanced. Quite naturally, the case of unbalanced 
population was first addressed a few years only after the BCS proposal. The question 
is of fundamental interest because the spin polarization opposes to the BCS pairing 
and eventually destroys superconductivity. Clogston |Tl and Chandrasekhar [12] 
independently considered the case where the spin polarization originates from the 
Zeeman coupling of the electron spin to an external magnetic field B. They had 
in mind a situation where the Meissner effect is absent: the field penetrates the 
complete metallic sample and the orbital coupling to the vector potential is negligible. 
This paramagnetic limit is characterized by the chemical potential difference ft = 
(/i^ — /i|)/2 = jB/2 between the two spin species. Here 7 denotes the gyromagnetic 
factor for conduction electrons. Clogston and Chandrasekhar found a first order 
transition at /i = Ao/^/2 above which superconductivity disappears and the normal 
state takes over the BCS ground state. Aq is the zero temperature gap for /i = 0. This 
result is simply obtained from the comparison of Gibbs energies, E = Eo — NoAq/2 for 
the unpolarized BCS state ([1]), and E = Eq — NqJI^ for the partially polarized normal 
state. ^0 is the energy of the normal state in the balanced case and A^o the density of 
states at the Fermi level. In fact, protected by the energy gap, the BCS state ([T]) is 
not modified when ft ^ 0. In particular, spin populations remain equal and the order 
parameter A remains equal to Aq. The first order nature of the transition is thus a 
consequence of the inability of the BCS state to react to a polarization constraint. The 
energy for flipping a spin down electron varies like Aq — ft with the chemical potential 
difference /i, but remains nevertheless positive at the transition. 

Unfortunately, for superconductors the Clogston- Chandrasekhar (CC), or PauH, 
limit corresponds in practice to very high magnetic fields. In most metallic compounds, 
superconductivity is limited by orbital pair-breaking effects and the corresponding 
upper critical field is much smaller than the Pauli limit. Orbital effects also lead to 
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Figure 1. Pairing in an imbalanced Fermi gas. a) In the classical BCS scenario, 
particles at the Fermi surface form zero center of mass momentum Cooper pairs. 
When the density of spin up and spin down particles are different, the two Fermi 
surfaces are mismatched and this mechanism must be modified, (b) FFLO: If 
particles with opposite spins form pairs at the surface of their respective Fermi 
seas, the Cooper pairs acquire a momentum, which can be interpreted as a spatial 
modulation of the order parameter, (c) Sarma: forming zero momentum pairs 
implies opening a gap inside the majority Fermi sea. 

a different physics with the emergence of a vortex lattice below the upper critical 
field. Pauli limited superconductors are rare and therefore correspond to specific 
conditions where the orbital coupling is suppressed, for instance 2D layers with an 
in-plane field. For more detail, we refer the reader to the concise but exhaustive 
review where the competition between Zeeman and orbital coupling is discussed in 
type-II superconductors. 

The details of the CC first order transition were later investigated by Sarma [13]. 
For ft > Ao/2, the normal state is locally energetically stable and an additional 
unstable solution of the gap equation appears. This intermediate unstable phase, often 
called the Sarma (or breached-pair [TT, Fig. [Hb) phase in the literature, connects the 
normal state when /i = Ao/2 to the fully paired BCS state A = Aq when ft = Aq. In 
the parameter range Ao/2 < ft < Aq, the variational energy as a function of A exhibits 
the typical features of a first order transition in the framework of Landau theory: two 
local minima located at A = and Aq, i.e. the normal and BCS phases, separated by 
a local maximum (the Sarma phase) at A^ = Ao(2/i — Aq). The Sarma phase shows 
interesting properties [151 IE] that have raised some interest ^171 [THl [191 EHl El] . Since 
A < ft in this phase, a significant number of Cooper pairs can be broken and the 
system becomes polarized in contrast with the standard BCS state [22]. Hence the 
Sarma phase possesses at the same time a superfluid component and gapless fermionic 
excitations between two Fermi surfaces. It is nevertheless an unstable phase [T3l [T7] . 
It has been noted [14 that the Sarma phase is a local energy minimum when the spin 
populations, and not the chemical potentials, are imposed. However the phase is then 
unstable against phase separation [23l [24] . 

Of course, these conclusions strictly hold only in the weak coupling limit. At 
unitarity, where the scattering length is infinite. Quantum Monte Carlo calculations 
have shown [25| that the gapless (Sarma) phase and the phase-separated BCS-normal 
state mixture are nearly degenerate in energy. The Sarma gapless phase was also 
found [20] to be stabilized in the presence of an optical lattice at intermediate coupling. 
In any case, no experimental evidence of a Sarma-type phase has been found so far. 

An original idea was put forward independently by Fulde and Ferrell [26] , Larkin 
and Ovchinnikov [27]: the lack of polarization flexibility of the BCS wavefunction 
can be improved by giving a finite momentum to the Cooper pairs (Fig. [He). The 
resulting order parameter A(r) is spatially varying. The finite momentum moves the 
Fermi surfaces of the two spin species with respect to each other, bringing them closer 
in one direction. This effect improves the pairing but costs some kinetic energy. The 
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displacement of the Fermi surface also leads to pair breaking and allows the system 
to polarize, thus relaxing part of the polarization constraint. As expected, the pair 
breaking is enhanced in the regions where the order parameter A(r) vanishes. There 
is often some confusion in the literature on the definition of the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) phases. They do not simply reduce to an helicoidal structure, 
A(r) = Ae*^ '^, or even a sinusoidal one, for instance A(r) = Acosq • r. They 
correspond generally to an inhomogeneous order parameter with a periodic structure 
that does not contain vortices. The characteristic length for the spatial variations 
is the coherence length of the superconductor. The phase transition between the 
normal state and the superfiuid FFLO phase was calculated in Refs. [28l[29^ 30j using 
the quasiclassical Eilenberger equations jSH [32l [33]. This is a first order transition, 
located at /i = O.TSIAq (the corresponding spinodal line is for ft = 0.754Ao), and 
it slightly extends the superfiuid domain when compared to the homogeneous BCS 
state as shown in Fig. [8l The order parameter has a cubic structure at the transition 
and evolves into a square and then a one-dimensional pattern as the temperarure is 
increased [34 . 

Interestingly, the structure of the FFLO phases as one penetrates inside the 
superfiuid remains an open problem, even at weak coupling. More is known on the 
transition between the homogeneous standard BCS phase and the FFLO phase. In 
analogy with magnetic systems, the superfiuid can form domain walls corresponding 
to a sign change of the order parameter over a spatial size on the order of the BCS 
coherence length. Far from the wall, the order parameter is constant. In addition 
to the usual gapped single-particle excitations, the wall shelters states with energies 
below the gap [35l [36] . Filling these states with spin up particles releases part of the 
polarization constraint, while the bending of the order parameter within the domain 
wall costs some condensation energy. A second order BCS-FFLO transition j37 l[38l [39] 
occurs when the wall energy becomes negative: a one-dimensional array of walls is 
formed, starting with infinite periodicity, which alternatively reverses the sign of the 
order parameter. The transition is located ^ dit fi 0. 666 Aq. As the chemical 
potential difference fi increases, the periodicity decreases merging the order parameter 
shape to a sinusoidal form, typical of FFLO phases. An analytical solution to this 
domain wall structure was obtained in the one-dimensional case [40l[4Tl[42] . Additional 
information on the physics of FFLO phases can be found in the complementary 
reviews [3 [43] [44]. 

In the context of cold atoms, the FFLO phases have been extensively studied 
theoretically along the BEC-BCS crossover and especially close to unitarity These 
various studies include mean-field calculations [H] [46l [47l [48] [49] , large N approach [50] 
or Density Functional Theory (DFT) [5T]. Including the trapping potential beyond 
the local density approximation has been argued to be necessary for describing the 
FFLO spatial oscillations [52l [53l [l^ . Nevertheless the same oscillations in the order 
parameter were interpreted in Ref. |54| to be a finite-size effect. With or without 
the trap, none of these approaches is completely exact and the results are therefore 
expected to be qualitative. Most works predict that the FFLO phases exist outside 
the weak coupling domain but there is no consensus on the value of the interaction 
strength for which they disappear. 
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3. Polarized Fermi gases 

3.1. Fermionic superfluidity with ultra-cold gases 

As mentioned earlier, Meissner effect prevents reaching the Pauh hmit in most 
superconductors. As a consequence, the Clogston-Chandraskhar instabihty scenario 
or the existence of Fulde-Ferreh-Larkin-Ovchinnikov phases were not unambiguously 
tested experimentally before the observation of superfluidity in ultra cold Fermi gases. 
These experiments were made possible by the development of techniques of laser 
cooling and trapping of atoms that led to the observations of the first Bose-Einstein 
condensates of alkali vapors in 1995 [55l [56]. Building on this breakthrough, a new 
generation of experiments was initiated soon after to cool down fermionic atoms and 
quantum degenerate Fermi gases were observed for the first time in 1999 by the group 
of JILA [57 . Fermionic superfluidity was obtained in 2003 [58l [59l [60l |6T] and allowed 
for the exploration of BEC-BCS crossover physics, a theoretical scenario bridging the 
gap between the Bardeen- Cooper- Schrieffer mean- field theory describing the behavior 
of weakly attracting fermions (scattering length a small and negative), and the strongly 
attractive regime (a small and positive) where the system behaves as a Bose-Einstein 
condensate of tightly bound dimers. Remarkably, although the theoretical foundations 
of this crossover physics had been laid in the early eighties by the pioneering works 
of Nozieres, Schmitt-Rink and Leggett [31 H], its first experimental confirmation 
was only made possible by the possibility of tuning interactions in ultra-cold atom 
vapors using Fano-Feshbach resonances. Another asset of cold atoms against classical 
condensed matter systems is the long spin relaxation time which offers the possibility 
of controlling spin populations using radio-frequency fields or optical pumping and 
keeping spin imbalances for long times. In addition to Feshbach resonances, several 
unique investigations tools were also developed and in the following we briefly review 
the main progresses that were achieved in the first years following the observation of 
fermionic superfluidity (note that we focus here on the tools used latter on for the study 
of spin imbalanced gases and as a consequence, we leave out important experiments 
such as the pair projection method [621 |63j or the study of collective modes [64] and 
let the reader refer to the more exhaustive review [1 ). 

3.1.1. Scattering length and Feshbach resonances Although ultra-cold gases are 
highly dilute, interactions still play a crucial role in their properties. Nevertheless, the 
description of interactions can be simplified by Pauli principle and the low temperature 
of these systems. Indeed, at ~ 1 /iK, the thermal wavelength of an atomic vapor is 
^0.1 /iK, and is therefore much larger than the typical range of interatomic potentials 
(~ 1 nm). This means that matter waves do not resolve the exact details of the 
potential, which can therefore be described as a zero-range contact potential. In 3D, a 
Dirac potential is highly singular and leads to divergences which need to be regularized. 
Several strategies can be followed, for instance the use of a non-local pseudo potential 
defined by 

Vpseudo{r) = gdr{r-)5{r), 

where the coupling constant is usually written as ^ = Anh^a/m and a is the scattering 
length ^. In the case of a two-body problem, this pseudo-potential is characterized 
by the following properties: 
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Figure 2. Left: At low energy, the short range interatomic potential is 
characterized by the scattering length a. A weakly bound state of energy 
El = —h'^/ma'^ is stable for a > 0. |a| = oo is associated with a universal 
regime known as the unitary limit. Right: the scattering length can be tuned 
using a magnetic field (here scattering length of ^Li between the two spin states 
|F = 1/2, mp = ±1/2} of the ground state hyperfine manifold). The divergences 
are known as Feshbach resonances and are associated with the resonance between 
a weakly bound state and scattering states. Fermions being stable close to a 
Feshbach resonance, it is possible to study the crossover between strongly and 
weakly attractive regimes. 

The scattering amplitude for a relative momentum k is given by f{k) = —a/ (1 + 
ika). We see that the interpretation of a is here straightforward, since it corresponds 
to the low energy Hmit of the scattering amplitude, and is thus associated with the 
scattering cross-section. An important point is the limit \a\ = oo, where the scattering 
length diverges, but the scattering amplitude stays finite and reaches an universal value 
1/ik. It can be demonstrated that this maximum value is imposed by the unitarity 
of the ^-matrix describing scattering events, and this regime is therefore called the 
unitary limit. 

For a > 0, the two body-potential possess a bound state of energy = 
—h^ /mo?. We see here that the unitary limit is associated with the disappearance 
of the two-body bound state, a generic feature which can be interpreted in terms 
of scattering resonances between incoming free particles and low-lying bound states. 
This connexion between the sign of the scattering length and the properties of weakly 
bound states suggests that the regime where a is small and positive corresponds 
to a strongly attractive potential, while on the contrary, the a negative and small 
regionn corresponds to a weakly attractive regime where the interatomic potential is 
too shallow to overcome quantum fiuctuations and maintain a bound state. 

In cold atoms, the scattering length can be be tuned by imposing an external 
magnetic field tuning the position of the bound state, the so-called Feshbach resonance 
j65l l66[ WJ\ [68j . Contrarily to bosonic system which are unstable when a diverges 
[69l[70], it could be demonstrated both theoretically [71 and experimentally [72 |l73[ [74 ] 
that Pauli principle was actually suppressing inelastic collisions in the regime of large 
scattering lengths. This unique property allows one to probe on a single atomic system 
the full crossover between strongly and weakly attractive domains (Fig. [2]). 
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3.1.2. Density profile and equation of state. Due to the external optical or magnetic 
confinement, the density profile of an ultra-cold gas is inhomogeneous and complicates 
comparison with theories usually developed for homogeneous systems. However, in 
most cases, the size of the system is large enough to define a mesoscopic length 
scale where the properties of the cloud can be described locally by the equation 
of state of an homogeneous gas. This is the Local Density Approximation (LDA). 
The local chemical potential /icr(^) can be defined, as the energy necessary to 
add a particle of spin cj at position r. The chemical potential /ihom[^t' ^i] 
homogeneous gas with the density Ua- is added to the trapping potential Ua- such 
that fiaif') = fJ^hom[n^ (r) ^ ni{r)] + Ua{r). At equilibrium, the local chemical potential 
jiia{r) is homogeneous and the density profile can be calculated from the resolution 
of the equation /ihom[^t(^)' ^^(^)] + ^cr(^) = where /io is a constant and is the 
global chemical potential of the system. 

Let us consider the special case of a spin balanced fermionic superfluid at zero 
temperature, with /j,^ = /j,^ = /j,. Dimensional analysis shows that the equation of 
state /ihom(^) takes the form 

/ihom(n) = Erfil/kra), (2) 

where = (67r^n)^/^ is the Fermi wave- vector of an ideal spin-polarized Fermi gas of 
density n, Ep = h'^kp/2m is the associated Fermi energy and / is some dimensionless 
function that must be determined from a microscopic theory and interpolates between 
the BEC and BCS sectors of the crossover. At unitarity (|a| = 00), the macroscopic 
equation ^ simplifies greatly since we have 

/^hom = £,Ef^ (3) 

where ^ = /(O) is just some numerical factor [75 . Remarkably, we observe that, 
up to this numerical factor, this is exactly the equation of state of an ideal gas: 
in other words, although we are considering the unitary regime, where interactions 
induce strong quantum correlations at the microscopic level, the macroscopic static 
properties are as simple as that of an ideal gas. In particular, in a harmonic trap 
the density profile is readily obtained using LDA, and the radius of the cloud is 
simply R = i^NiFGv^, where i^NiFG is the radius of a non interacting Fermi gas 
with the same atom number and in the same trapping potential. This relationship 
suggests a straightforward way to determine the interaction parameter ^, by simply 
measuring the ratio R/RmFG- This very scheme and some variations have been 
applied experimentally in situ or after time of flight and theory [76l [771 1211 EHl [79] 
and experiments [8(3 [HU [821 [83] now converge towards the value ^ = 0.42(1). 

However, the simple scaling law presented above is valid only at unitarity and 
comparison between theory and experiments in the crossover is more involved. Indeed, 
atomic density profiles are obtained in practice by absorption imaging which gives 
access to an integrated two-dimensional density n2D{x^ y) = J dz n(x, y, z), where the 
hue of sight is taken along z. To obtain insight on the local properties of the cloud, it is 
thus necessary to deconvolve this integration. In the case of a cylindrically symmetric 
trap, this can be achieved by using a data processing technique inspired from medical 
imaging technology based on a mathematical operation called Abel transform [84J and 
applied recently in the case of imbalanced Fermi gases [85 . An alternate scheme to 
measure the grand canonical equation of state of the homogeneous gas was proposed 
recently [86 . Indeed, let us consider the double integrated density of spin i atoms 
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defined by ni{z) = J dxdy n{x^y^ z). Using Gibbs-Duhem relation = d^^P relating 
pressure P and density n^, we can write that 



Due to the spatial dependence of the chemical potentials, LDA allows one to turn 
the spatial integral into an integral over (ii reading 



where uo'^ = uJx^y is the transverse trapping frequency. Integration is straightforward 
and yields 



Eq. implies that one gets access to the equation of state P(/ii,T) of the 
homogeneous gas simply by measuring the doubly integrated density profile in the 
trap. This relation is very general and can be used in a harmonic trap as soon as local- 
density approximation is valid, for instance in the presence of a spin imbalance or in an 
optical lattice. Remarkably, the equation of state P(/ii, T) contains all the macroscopic 
information on the system, since using Gibbs-Duhem relation, it is possible to extract 
missing thermodynamic quantities, such as atom density or entropy. This scheme as 
been applied for the first time in [S'f to the case of the finite temperature equation 
of state of a balanced Fermi gas, but can also be used to probe the properties of 
imbalanced systems, as we will see later. 

3.1.3. Radio- frequency spectroscopy. One of the main features of fermionic 
superfluidity is the presence of a pairing gap in the excitation spectrum that constitutes 
the most direct manifestation of Cooper pairing responsible for the onset of quantum 
order in these systems. In cold atoms, this quantity can be accessed by spectroscopic 
tools: for instance, it is possible excite atomic spin degrees of freedom using radio- 
frequency field, and measure the shift of the resonance induced by interatomic 
interactions. Several groups have measured the binding energy of fermion pairs in 
ultra-cold Fermi systems using this scheme [88] , starting from pioneering experiments 
by the groups of JIL A and MIT on quantum degenerate Fermi gases [89^, [90] , to the 
recent measurement of the spectral function of a fermionic superfluid at JIL A |9T| . 
Nevertheless, the interpretation of experimental data is non-trivial since interactions 
between the superfluid and the final spin state can modify strongly the response of 
the system. On the theory side, intense effort has been devoted to the understanding 
of the role of interaction in rf spectroscopy |50l [93l [94l [95], and in experiments, it 
has been pointed out that initial and final spin states could be chosen so as to reduce 
these spurious effects and permitted an unambiguous measurement of the excitation 
gap of the system as demonstrated in [96| . 

3.1.4- Vortices and superfluidity. Although in BCS theory pairing and superfluidity 
arise at the same time, this property is not general. In most cases the two phenomena 
are decoupled: for instance, in the case of a Bose-Einstein condensate of deeply bound 
molecules, the formation of the dimers takes place at a temperature much higher 




(4) 
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than the condensation critical temperature. In practice, the most convincing probe of 
superfluidity in ultra-cold Fermi gases was obtained by the group of MIT by stirring the 
cloud using a far detuned laser beam and observing the formation a triangular array 
of quantized vortices that can be considered as a the smoking gun for the existence of 
a complex order parameter characterizing the superfluid phase [97J. 

3.2. Ultra- cold Fermi gases with imbalanced spin populations 

Once balanced fermionic superfluidity is achieved with ultracold systems, the 
experimental study of spin polarized Fermi systems is straightforward. The groups 
of Rice University and MIT prepared simultaneously for the first time ultra-cold 
Fermi gases with imbalanced spin populations [82^^^, and, despite some discrepancies 
between their results, the two experiments support Clogston-Chandrasekhar's original 
proposal that the superfluid phase is robust against spin polarization. Indeed, both 
groups observed that, in the presence of an imbalance in the spin populations, the 
cloud phase separates along a shell structure. Using Abel reconstruction they observe 
that a fully paired core with equal spin densities is surrounded by a "magnetized" rim 
with imbalanced spin populations [99, 100 . By rotating the cloud, it was in addition 
observed that only this central core could sustain quantized vortices, indicating that 
the central phase is superfluid and that the magnetized phase is not [98 . This 
very simple observation constitutes a crystal clear confirmation of the Clogston- 
Chandrasekhar hypothesis. Indeed, for cold atoms, a polarizing magnetic field is 
equivalent to a chemical potential mismatch between the two spin species. In the 
case where the trapping potential is the same for the two spin states (which is the 
situation realized experimentally), LDA shows that the chemical potential difference 
ft = (/i^ — /i^)/2 is constant in the cloud. Since the gap A decreases with density, the 
CC scenario predicts that an unpolarized superfluid should occupy the center of the 
trap, as long as A stays roughly larger than /i, which is qualitatively in agreement 
with experimental findings. 

Despite the important similarities in the structure of the central core obtained 
by the Rice and MIT experiments, major discrepancies exist on the structure of the 
normal phase at the rim of the trap. Indeed, while Rice's group observes a single fully 
polarized phase where only majority atoms are present, MIT obtains a richer phase 
diagram, with a normal component divided between a fully polarized region similar 
to the one observed at Rice, and a mixed phase where the two spin component are 
present (Fig. [3]). A first consequence of this difference is that apart from the extreme 
polarization case = 0, a superfluid core is always present in Rice's observations. 
This particular feature contradicts MIT measurements, which observes that above a 
polarization Pc = (A^^ — Ni)/{N^ + N^) (Pc ~ 0.75 at unitarity, confirmed by recent 
ENS measurements), the core is no longer superfluid, and only the two normal phases 
are observed. 

4. Zero temperature phase diagram 

In this section we present theoretical interpretations of the experimental results 
outlined above. Since most experiments were performed at unitarity, we will focus 
first on this regime, and then we will later provide an extension to the whole crossover. 
A key question we wish to answer is whether a stable partially polarized normal phase 
exists. This issue constitutes the main difference between the experiments of Rice and 
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Figure 3. Density profile of a unitary spin imbalanced Fermi gas. From top to 
bottom: n-^ (majority, red), fid = n-f — (green) and (minority, blue). At 
low temperature, the density profile presents a shell structure, with a fully paired 
superfiuid at center, surrounded by a partially polarized normal phase, and at 
the outside a fully polarized ideal Fermi gas of majority particles. Note that the 
fiat profile of the integrated density difference Ud is the signature of a fully paired 
core [ToT] . 



MIT and we shall show that this question can be answered by the study of a simpler 
problem, namely the calculation of the energy for an impurity immersed in a Fermi 
sea. Although rather academic at first glance, the resolution of this problem will 
constitute the backbone for the description of the normal phase in MIT's experiment. 



4.1. Scaling law and unitary phase diagram. 

Due to the presence of a trap in experiments, it is more convenient to describe the 
behavior of the cloud in the grand canonical ensemble, where the thermodynamic 
properties of the system are fully characterized by the equation of state P(/i^^^^|,T) 
|lQ2j . Using the same dimensional argument as for the balanced unitary gas, this 
equation of state at unitarity and T = can be recast as a function of a single 
dimensionless function h such that 

where Pq is the pressure of a single component ideal gas [1Q3| 1104] . 

A sketch of h corresponding to Rice's observations is displayed in Fig. IHa. On 
the one hand, the equation of state in the normal phase is that of an ideal gas of spin 
up particles, hence h = 1. On the other hand, the knowledge of the equation of state 
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Figure 4. Sketch of the equation of state of an imbalanced Fermi gas. From 
the thermodynamics relation Q = —PV, the stable phase is that with the highest 
h. (a) Rice experiment. Only two phases were observed: a fully polarized ideal 
gaz (black horizontal line) and a fully paired superfluid (red line). They cross at 
?7o ~ —0.1 which indicates the position of the Clogston-Chandraseckar limit in this 
situation, (b) MIT experiment. The dashed and dot-dashed line give a sketch of 
stable and instable intermediate normal phases. Here r]c is the normal-superfluid 
threshold, and A the boundary between the partially and fully polarized phases. 
Graphically, we observe that the stability criterion for the intermediate phase is 
A < ?7o. Moreover, A/j.^ can be interpreted physically as the energy of the last 
minority atom removed from the majority Fermi sea, and is thus closely related 
to the Fermi-polaron problem. 



of a balanced superfluid, Eq. (j3j), yields in the superfluid core h = (2^)~^/^(l + r^)^/^ 
jlQ5j . In this diagram, the CC limit corresponds to r^o = (2^^^^ — 1 ^s.^ —0.1 where 
the two curves meet. Indeed, since pressure is related to the grand potential Q and 
the volume V by the relation Q = —PV^ the most stable phase at given chemical 
potentials is the one with the highest pressure. For rj > rjo the most stable phase is 
the fully paired superfluid, while the fully polarized ideal gas is favored for r] < rjo- 

Let us now consider what would be the phase diagram corresponding to MIT's 
observations, as shown in Fig. [Hb. The intermediate mixed phase is a partially 
polarized normal phase that intersects the fully polarized Fermi gas at 77 = A. 
Physically, the corresponding phase transition is interpreted as the removal of the 
last spin down particle from the partially polarized normal phase. Therefore the 
dimensionless parameter A is determined from the energy £^ of a single minority 
atom immersed in a Fermi sea of majority particles. By definition, this energy is the 
chemical potential /i^, while ja^ is the Fermi energy Ef of the spin up Fermi sea, 
hence A = E/Ep |1Q4[ 1103] llQ6j . The value of A is crucial to determine the relevant 
scenario for the density profiles. For A> rf^^ the mixed phase is unstable (dot-dashed 
line in Fig. [Hb), always preempted by the superfluid phase, and RICE's features are 
recovered. However, for A < 7^0, the mixed phase (dashed line) extends between r] = A 
and T] = r]c^ in agreement with MIT's observations. 

4.2. The Fermi polaron. 

In the previous section, we have seen that the stability of a partially polarized phase 
could be related to the physics of an impurity immersed in a Fermi sea, the so- 
called Fermi-polaron named after the solid state physics polaron that describes the 
interaction of an electron with a bath of (bosonic) phonons. Despite its simplicity (one 
impurity immersed in a non interacting Fermi sea), the lack of a small parameter in 
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this problem makes it strongly interacting and essentially non-pert urbative. Monte- 
Carlo methods are a priori required to calculate the energy spectrum of the impurity 

However, to demonstrate the stability of the normal phase observed at MIT, we 
just need an upper bound to this energy. A variational method is thus sufficient 
to settle the issue at stake. In the polaron picture, an impurity immersed in some 
medium forms a quasi-particle composed of the bare impurity dressed by excitations 
of the surrounding environment. In the present case, elementary excitations of a Fermi 
sea are particle-hole pairs that shroud the spin down atom to form the Fermi-polaron. 
Since the impurity is alone, we can assume that the effect of the impurity on the Fermi 
sea remains weak, and we can try to work in a subspace where a single particle-hole 
pair is created by the presence of the spin down atom. Variational equations can be 
worked out analytically, leading for small p to a dispersion relation 

Ek = AEf + 7— + 
2m* 

analog to that of a free particle with renormalized parameters A ~ —0.6 and 
m* = 1.17m for a = oo jlQ3[ 1108] . Let us stress that this result is variational, 
and gives us access to an upper bound of the energy. The true value of A is then 
strictly speaking smaller that -0.6. Nevertheless, this bound is sufficient to conclude 
that A is indeed smaller than rjQ and to demonstrates rigorously the existence of a 
partially polarized phase between the fully polarized and fully paired sectors of the 
phase diagram. 

Even more remarkably, it appears that the variational values found here are 
strikingly close to the actual predictions A = —0.59(1) and m* = 1.09m for the 
variational Fixed Node Monte- Carlo |106j and A = —0.61(1) and m* = 1.20 (l)m for 
diagrammatic Monte-Carlo [107 . This coincidence can be explained by the relatively 
weak probability of excitation of a particle-hole pair. Indeed, in the variational 
calculation this probability is only ~ 25%, despite the strength of interactions at 
unitarity . This is also conffimed by the moderate modification of the polaron mass 
which is close to its bare value, and also by a systematic expansion of the polaron 
energy as a function of the number of particle hole-pairs excited which demonstrated 
that the associated series in kp was converging quite fast [107' "109". 

Let us now discuss the equation of state of this intermediate normal phase. As 
ffist proposed in |106j . if one assumes that the impurities keep their fermionic nature, a 
natural description of their collective behavior at low density of spin down particles is 
that of an ideal Fermi gas of polarons, a picture conffimed by Monte Carlo simulations 
of the equation of state. In this regime, it was shown in Ref. j87 |lll0j that the pressure 
of the normal mixture is the sum of the pressures of an ideal gas of majority atoms 
and an ideal gas of polarons, i.e. 

Despite its seemingly non-interacting form, Eq. ([7]) does include interactions between 
polarons as a result of the fi^ dependence in the second term. Translating Eq. ([7]) into 
the canonical ensemble, an exact relation between the dominant polaronic interaction 
and the parameter A was derived in Ref. [110 , in relatively good agreement with fixed 
node Monte-Carlo simulations (within a few percents). 
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We observe that the equation of state of the normal phase meets that of the 
superfiuid at a value r]c ~ 0.065. It corresponds to the maximum chemical potential 
mismatch that the superfiuid can sustain before turning into the normal state. 
This determination of the Clogston-Chandreskhar limit at unitarity allows one to 
reconstruct the density profile in a trap, as well as calculate the critical polarization 
at which the superfiuid core vanishes in MIT's experiments [111", yielding Pc ~ 0.77, 
in close agreement with the value measured at MIT and ENS. 

4.3. Extension to the crossover and the polaron-molecule transition. 

In the previous section, we have shown that the normal phase at unitarity could 
be described quantitatively by a Fermi liquid theory involving Fermi polarons whose 
physical properties could be obtained with a rather good accuracy from a simple 
variational calculation. This picture is still valid in the BCS sector of the crossover, 
since by construction, the variational ansatz used to describe the Fermi-polaron 
coincides with second order perturbation theory in the limit a ^ 0~ . In the BEC 
sector, however, the Fermi liquid scenario must break down at some point since in the 
far BEC limit, we expect the impurity to form a tightly bound dimer with one majority 
atom, giving rise to a bosonic, rather than fermionic, behavior of the dressed quasi- 
particle. This molecule/polaron transition was unveiled by Monte-Carlo simulation of 
the impurity problem [107 . This calculation shows that the transition takes place at 
l/kr^ci = 0.91(2), and that above this threshold the energy of the polaron is accurately 
described by that of a point-like boson interacting with the Fermi sea by a mean-field 
energy ^ad^t' where the atom-dimer coupling constant ^ad is characterized by the well 
known atom-dimer scattering length aad ^ 1.18a |112j . This unexpected robustness 
in the validity of the mean- field approximation was explained later on by an extension 
of the variational description of the polaron to the molecular sector |113[ I114[ I115j . 
Similarly to the BCS/unitary case, the quasi-particle is described as a dimer dressed 
by single particle- hole pair of the majority Fermi sea. Due to the composite nature of 
the dimer, the calculations are a little more involved than in the fermionic case, but 
can be reduced to a single integral equation that is solved numerically. The energy 
displayed in Fig. [5la demonstrates a remarkably good agreement between Monte-Carlo 
simulations, variational calculations and rf experiments discussed below. 

From bosonic to fermionic along the crossover, the change in the nature of the 
impurity immersed in the Fermi sea of majority atoms has strong consequences on 
the structure of the phase diagram. Indeed, while no qualitative difference is expected 
between unitarity and BCS sectors (essentially the same three phases that we discussed 
above), the BEC regions is dominated by a completely different physics \116\ . As 
suggested by the molecular scenario, we can describe the system in the BEC limit as 
a mixture of a Bose-Einstein condensate of dimers and a Fermi sea of excess fermions 
using the mean-field equation of state |117| 1118] 

E f N ^ ^ fa 2\2/3 5/3 , ^dd^^^ , 

yi^F.nB) = -^nB^-—{67r) + + ^ad^B^F + ... 

with ub = \^ the number of bosonic dimers, np = — is the number of 
free excess fermions and ^dd is the coupling constant describing s-wave dimer-dimer 
interactions and characterized by the scattering length add = 0.6a (TTJ I119j . In 
contrast with the fermionic physics at play in the BCS region of the phase diagram 
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where the gap maintains a perfect pairing, for n^? 7^ 0, this equation of state describes 
a polarized superfluid. For weak enough repulsive interactions between bosons and 
fermions, the Bose- Fermi mixture is stable, while strong interaction close to the 
Feshbach resonance drives a first order phase separation [ 1171 1116^ I120[ I121j between 
a mixed Bose- Fermi phase and a fully polarized Fermi gas. The Fermi gas is made of 
majority fermions repelled from the Bose-Fermi mixture. A linear stability analysis 
yields the stability criterion 



that is obtained from the positivity of the compressibility matrix d^.^.E. Using the 
values of aad and add given above, we finally get that an atom/dimer phase separation 
takes place for l/kp^a < 1.7 when ^ (For comparison, the mean field prediction 
is l/kp^ci 1.88 |122| ). This onset of a first order transition gives rise to a tricritical 
point where the order of the transition changes and that will be studied in more detail 
in the paragraph dedicated to finite temperature phenomena, see below. The linear 
stability analysis also predicts a critical end-point at zero polarization after which the 
bosonic superfiuid can no longer sustain a finite polarization. All additional majority 
fermions are therefore expelled. The corresponding stability criterion is given by 



4) [2m) ^''"^ 



-t = nB<(^) {^] (6.Yt^, (9) 



for = n^. The critical end-point is thus located at l/kp^a ~ 0.66 in fairly good 
agreement with the value l/kp^a — 0.53 obtained with Fixed Node Monte-Carlo |116j . 

The complete canonical phase diagram in the BEC-BCS crossover was first 
established from a Nozieres-Shmitt-Rink analysis revealing its most important 
features [122 . A more quantitative description was later provided by Fixed Node 
Monte-Carlo simulations [116 and in Fig. Ob, we present a grand canonical 
description of the phase diagram based on this work. 

From dimensional analysis, the structure of the phase diagram can be 
encapsulated in two dimensionless numbers: rj = /i^Z/i^ which characterizes the spin 
imbalance and S = h/ ^^2miM^a which for an ideal gas would be equal to l/kpa and 
thus measures the strength of interactions. Using Local density approximation. Fig. 
Ob can be used to predict the shape of the density profile in a trap. From equation 
/i^(r) = fi^ — U(r) and the definition of (5, it is indeed possible to express U{r) as a 
function of S{r). Substituting in 7^ = fi^/f^i, we then see that in a trap, the system 
follows a line rj = 1 — 4J^/i^ma^/fi^, with ft^ = (/i^ — /iJ)/2. Interestingly, this equation 
does not depend on the actual shape of the trapping potential which only plays a role 
in the relationship between chemical potentials and atom numbers. 

While the BCS side of the resonance has been thoroughly studied experimentally, 
the predicted features of the phase diagram on the BEC side have not yet been 
explored and await an experimental confirmation. The increase in the atomic loss 
rate towards the BEC region limits current experiments close to the unitary and 
BCS regions. In fact, even on the theoretical side, a clear understanding of the 
bosonic/fermionic crossover at large polarization has not been achieved yet. This 
is related to the question of whether the polaron/molecule transition is affected, or 
even preempted, by the demixion mechanism described above, which leads to the 
coexistense of a bosonic superfiuid and a single-component Fermi gas. Two scenarios 
can be envisioned. The first one is supported by Quantum Monte-Carlo calculations. A 
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Figure 5. Imbalanced fermi gas in the BEC-BCS crossover, (a) energy of 
an impurity immersed in a Fermi sea and molecule/polaron transition. Solid 
(dashed) hne, variational calculation of the polaron (molecule) [1141 11131 I115| . 
Solid circles: experimental results from MIT [ 123j . solid (empty) diamonds: 
Diagrammatic Monte-Carlo calculation of molecule (polaron) [107 . (b) phase 
diagram in the crossover in the grand canonical ensemble. The gray region 
indicates the superfluid sector. PP stands for partially polarized. Solid (resp. 
dashed) lines correspond respectively to first and second order transitions. Note 
that beyond 6 = 1.7, the normal/superfiuid transition becomes second order. The 
dotted lines correspond to the "path" followed in a trap. 



first order transition occurs, at [116" l/kp^ci = 0.73 in the impurity limit 0, where 

the polaronic normal phase separates into a polarized paired superfluid immiscible 
with the remaining majority fermions. This transition hides the polaron/molecule 
transition, that is located [107 further at l/kp'^a = 0.91(2). The possible uncertainty 
in Monte-Carlo calculations and the proximity of these two transition points allow to 
imagine an alternative scenario where the two points indeed coincide. In that case, the 
polaron/molecule transition would be accessible and immediately followed by phase 
separation on the BEC side. Interestingly, even the nature of the polaron/molecule 
phase transition has not been unambiguously determined and a smooth transition 
remains possible. 

4.4- Experimental characterization of the Fermi polaron. 

The review of imbalanced Fermi gases presented in the previous sections has shown 
that the behavior of the partially polarized normal phase could be essentially 
understood from the properties of the Fermi polaron. We show here that some of the 
experimental investigation tools developed for balanced Fermi gases can be applied to 
imbalanced situations. These studies confirm quantitatively some of the theoretical 
predictions presented above, in particular the spectrum of a single polaron or the 
equation of state of the partially polarized normal phase. 

4-4- ^- Radiofrequency spectroscopy. The parameter A that characterizes the 
transition between a fully polarized and a partially polarized normal phase was 
measured in Ref. [123 using radio- frequency spectroscopy of a strongly imbalanced 
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gas. When the imbalance is large enough (more specifically larger than Pc ~ 0.75 at 
unitarity), no superfluid core is observed, and all minority atoms are located in the 
normal component. The rf spectroscopy of spin down atoms therefore gives access 
to the energy shift AEp^ induced by the presence of the majority Fermi sea. These 
results are presented in Fig. [5] and confirm remarkably well the theoretical approaches 
developed to describe the crossover. In addition, the measurement of the spectral 
weight of the radio- frequency resonances allowed for a characterization of the fermionic 
nature of the polaron and demonstrated that beyond l/kp^a ~ 0.74(4), the impurity 
was loosing its fermionic nature. This value is rather far from the molecule/polaron 
transition but is by contrast strikingly close to the prediction of the Fixed Node Monte- 
Carlo for the transition between the partially polarized superfluid and normal phases. 
It therefore suggests that, as discussed above, the single impurity molecule/polaron 
transition is preempted by the phase separation between the ideal Fermi gas and the 
polarized molecular superfluid. 

4.4-2. Dynamics and effective mass. To measure the effective mass of the polaron, 
it is necessary to study the dynamics of an imbalanced Fermi gas. As proposed in 
[106 , this can be done by the study of the collective modes of the mixture. Indeed, 
the energy of a trapped polaron is given by the equation 

2 2 

E{r,p) = AEF^ir)+£-^+U{r) = AEF^{0)+£-^+{l-A)U{r), (10) 

where we have used the LDA condition Ep^ir) = Ef^{0) — U{r). In the case of a 
harmonic trapping U{r) = y ^rnujfxf we see that the trapping potential of 

the polaron is still harmonic, with an effective frequency 

uj^/uji = ^(1 - A)m/m*. (11) 

Using the magnetic field dependence of the scattering length, it is possible to 
excite selectively the minority atoms by canceling the interatomic interactions. Such 
an experiment was performed in Ref. [124^ . The study of the axial breathing mode of 
an elongated imbalanced Fermi gas showed that at unitarity the oscillation frequency 
was about 1.17(2) times bigger than that of an ideal gas. From this measurement, 
Eq. (pTj) and the value A = —0.6 obtained from theory or rf measurement yield 
m* = 1.17(10)m, in agreement with theories of the polaron presented previously. 

4.4-3. Collective behavior. In the polaron picture, the partially polarized normal 
component is described by the equation of state of a mixture of two ideal Fermi gases 
(Eq. [5]). Using the general method described in section [3.1.21 it is possible to directly 
measure this grand-canonical equation of state by the analysis of absorption imaging 
pictures, as displayed in Fig. Oa [87 . A remarkable agreement is thus obtained 
between experiments and the Fermi liquid model elaborated for the properties of the 
Fermi-polaron quasi-particles. 

This procedure can be extended to the full crossover and, except in the far BEC 
region of the phase diagram where the polaron picture breaks down the equation of 
state ([7]) agrees with Monte-Carlo simulations |125) as well as experimental data [126]. 
This observation confirms that Fermi-polarons are indeed fermions whose effective 
parameters A and m* are well understood using the variational model described above. 



Ultra- cold Polarized Fermi Gases 



17 



2.2 



2.0 
1.8 
^ 1.6 
1.4 
1.2 
1.0 

-0.6 -0.4 -0.2 0.0 0.2 

Figure 6. Equation of state of a unitary imbalanced Fermi gas in the grand 
canonical ensemble. The mixed normal component spans 77 = [A, ijc] ~ 
[—0.6,0.065]. The dashed line corresponds to the polaron Fermi liquid Eq. [Jl 
while the dotted line is the prediction of the Monte-Carlo prediction [106 . Above 
r] = r]c, the equation of state is well described by that of unitary gas with 
= 0.42(1). Inset: density ratio n^/n^ calculated from the thermodynamic 
identity rii = d^-P. The critical value r] = rjc separates a superfluid with = 
and a partially polarized normal phase. Figure from [Sj. 

4-4-4- T^he bosonic sector The experimental exploration of bosonic region of the 
crossover is presented in jl27 j. This work first confirms the disappearance of 
the polaron phase for l/kpi > 0.75, in agreement with Monte-Carlo simulations. 
The analysis of the cloud density profile furthermore demonstrates that beyond 
this threshold the system behaves like the Bose- Fermi mixtures described above. 
Quantitative comparison between theory and experiment allowed for the measurement 
of the dimer/dimer and atom/dimer scattering lengths confirming the expected values 
ci^dd = 0-6a and aad = 1.2a. Interestingly, experimental data showed small deviations 
from the bosonic mean-field behavior that could be explained by the introduction of 
beyond- mean- field Lee-Huang- Yang term [128j . 

5. Finite Temperature phase-diagram 

The phase diagram discussed in the previous section displays intriguing features when 
extended to the finite temperature regime. Indeed, thermal fluctuations reduce the 
value of the gap A of the excitation spectrum and, according to the CC argument, 
weaken the robustness of the superfluid against spin polarization. We thus expect the 
normal-superfluid transition to follow a critical line r]c{Tc) in the {rj^kBT / ji^) plane, 
starting from 7?c(0) = 0.065 as discussed earlier, and ending at ri{Tc) = 1, for some 
critical temperature Tc for which the superfluid displays no longer any resistance to 
polarization. Since the phase transition takes place at zero imbalance, we identify 
Tc with the usual superfluid/normal second order transition that was was recently 
located at ksTc/ ja = 0.32 at unitarity jl29[ 11301 [87] . In the BCS weakly interacting 
regime, the position of this transition can be obtained from BCS theory. Remarkably, 
while at zero temperature the transition is first order, as discussed above, we expect it 
to become second order at zero imbalance. This suggests the existence of a trier it ical 
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Figure 7. (a): Sketch of the finite temperature phase diagram of the unitary 
imbalanced Fermi gas. The phase transition being first order at zero temperature 
(circle) and second order at zero polarization (diamond) implies the existence of 
a tricritical point located at (?7t,Tt), (triangle). Above this point, the critical line 
is second order (Dashed line), and first order below (Solid line), (b): Tricritical 
temperature in the BEC-BCS crossover. In the far BCS limit, Tt is proportional to 
Tc and is thus exponentially small. At l/kjp^a = 1.7, the tricritical temperature 
vanishes and coincides with the zero-temperature demixing threshold from Eq. [8] 
(diamond) (Data from [122| , courtesy of M. Parish). 



point somewhere in between as seen in Fig. ([7j) which outhnes the main features of 
the finite temperature phase diagram of a unitary Fermi gas. 

Interestingly, the Bose/Fermi behavior of the polarized Fermi gases on the 
BEC side of the resonance has similarities with ^He/^He mixtures |131j and 
the phase diagrams thus share common features. In particular, a first order 
phase transition between phases with different isotopic ratios were observed at 
low temperature in ^He/^He. Since this phenomenon forms the basis of dilution 
refrigerators, the physics of helium mixtures is of the highest practical interest and 
their thermodynamic properties around tricritical points have been studied extensively 
using Renormalization Group approaches in the early 70 's |132( I133| I134j . In this 
context, it was demonstrated that, within logarithmical corrections, the critical 
exponents in the vicinity of the second order /first order transition in 3D were given 
exactly by Landau theory. 

In the case of ultra-cold gases, renormalization group calculations have been 
used to derive fairly accurately the position of the tricritical point at unitarity |135j , 
and yielded Tt/Tp^ = 0.09. This value is in close agreement with experimental 
observations that measured its position by the disappearance of the density jump 
at unitarity associated with the first order nature of the normal-superfiuid transition 
a large imbalance [136 . A semi-quantitative theoretical extension to the BEC-BCS 
crossover has been proposed in [122J. Using Nozieres-Shmitt-Rink approximation, it 
is possible to show that the temperature of the tricritical point follows a curve which 
connects the unitary tricritical point to the zero-temperature demixing threshold (Eq. 
[8]) where, as mentioned earlier, the normal/superfiuid transition goes also from first 
to second order. 

The tricritical point bears an intriguing connexion with FFLO phases (Fig. [8|). 
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Figure 8. Mean field phase diagram including FFLO phases (BCS limit). This 
graph generalizes the phase diagram of Fig. \7\to the case where the FFLO phase 
is taken into account. The tricritical point corresponds also to the disappearance 
of the inhomogeneous superfiuid. 



Indeed, in the weak coupling BCS regime, the domain with FFLO phases shrinks as 
the temperature is increased from zero: the first order FFLO-normal state transition 
line and the second order homogeneous(BCS)-FFLO transition line meet exactly at the 
tricritical point P^l34j, located at Tt/Tc ^ 0.561, pt/^o — 0.608 for weak coupling, 
where Tc is the critical temperature in the balanced case. Above the tricritical 
point, the superfluid- normal transition is the standard second order transition to the 
homogeneous BCS phase. Interestingly, the Clogston-Chandrasekhar limit is always 
preempted by the FFLO transition and also ends at the tricritical point. This feature 
is specific to the weak coupling limit and other scenarios can be envisioned at 
stronger coupling [137] . 

6. Conclusion 

Experimental results obtained in the past years at MIT, Rice, and ENS have offered 
a simple picture of the phase diagram of an imbalanced Fermi gas with attractive 
interactions. Quite surprisingly, the main features can be captured with a remarkable 
precision by models based on simple concepts such as Fermi Liquid theory or polaron 
physics. Despite these undeniable successes, many open questions remain and will be 
explored in the near future. The most serious one is the discrepancy observed between 
Rice results and those of MIT and ENS. The absence of a normal mixture in Rice's 
observations is still an unsettled issue and several explanations have been proposed. 
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among which the role of trap anisotropy, atom number |10H 1138] or evaporation |139j . 
However, ENS experiments were performed in a regime close to that of Rice, and do 
shows features very similar to those of MIT. 

The second point is related to the existence of a FFLO phase, and more generally 
that of a the existence of a partially polarized superfluid. Although no systematic 
experimental study of this question has been performed in 3D, the remarkably good 
agreement between theory and experiment suggests that the FFLO phase can occupy 
only a small region of the phase diagram and/or do not influence much the macroscopic 
properties of the cloud. It was nevertheless suggested that, in the spirit of the 
celebrated Kohn-Luttinger theorem j l4Qj . at low temperature the polaron gas could 
form a p-wave superfluid |141j . but up to now, no experimental signature of this 
new phase was reported. It was known from earlier works |3Q1 HI] that a reduced 
dimensionality, in particular in ID, improves the Fermi surface nesting and therefore 
favors the emergence of FFLO-like features. In the last few years, a considerable 
theoretical effort was devoted to understand and predict the existence and fingerprints 
of the FFLO polarized superfiuid state in ID, combining bosonization approaches |142[ 
[143] , Bethe-Ansatz calculations of the Yang-Gaudin model [T44l [T45l [146], Monte- 
Carlo [Mil [UHl [Ha] and DMRG studies [Ml [HH [ISS [l53l MM MEi [US] of the 
attractive Hubbard model or of the Yang-Gaudin model, and a spin- density- functional 
theory [157^. Recently, very promising results where obtained in this direction by the 
group of Rice who trapped ultra-cold fermions in elongated quasi-lD tubes |158) . 
The study of the density profiles was in good agreement with theoretical predictions 
displaying FFLO-like behavior. 

Third, despite important progress in the understanding of static properties of 
imbalanced Fermi gases, very little is known on their dynamic behavior. One challenge 
is the understanding of the crossover between the hydrodynamic and collisionless 
regimes at low and high polarizations respectively. First experimental evidences 
were presented in |124j . but a unifying theoretical description is still missing, despite 
some attempts using Fermi-Liquid theory [ 159| 1160] . or the resolution of coupled 
Boltzmann/Euler equations [161 . 

Finally an interesting perspective is offered by the recent realization of quantum 
degenerate Lithium-Potassium Fermi mixtures |162H163] . In particular, recent results 
on the width of the Feshbach resonances |164j and the stability of the ^Li/^^K 
mixture |165j offer the opportunity to explore experimentally a new generalization of 
the BCS superfiuidity mechanism where mismatch between Fermi surfaces is created 
by mass imbalance, leading to a phase diagram even richer than the one presented 
here [ml [Ml |M1 [Ml [1691 [TTO] . 
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